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Abstract 

We perform measurements, numerical simulations, and quantitative comparisons with available 
theory on solitary wave propagation in a linear chain of beads without static preconstrain. By 
designing a nonintrusive force sensor to measure the impulse as it propagates along the chain, 
we study the solitary wave reflection at a wall. We show that the main features of solitary wave 
reflection depend on wall mechanical properties. Since previous studies on solitary waves have been 
performed at walls without these considerations, our experiment provides a more reliable tool to 
characterize solitary wave propagation. We find, for the first time, precise quantitative agreements. 
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Solitons are widely studied in physics because of their ubiquity in systems exhibiting non- 
linear propagation |l|. In a granular chain, theoretical and experimental evidence of solitons 
was first reported by Nesterenko ^ [1 U, 0, 0, Q| . Since Nesterenko's pioneering work, most 
of the experimental effort in the field has generally focused on the scaling laws for amplitude 
and speed of the solitons, 0,0]. It was recently reported 10] that identical and opposite 
propagating solitons do not preserve themselves upon collision and hence these are solitary 
waves rather than solitons. Several detailed numerical studies have been devoted to under- 
stand the interactions of solitary waves with a perfectly reflecting wall 0, H 0, Q 
and show that tiny secondary solitary waves are generated as a solitary wave is reflected 
off a wall jloL 0] . However, due to experimental difficulties, no close comparison between 
experiments and simulations has so far been established. Here inspired by Nesterenko's 
experiments 0, (J [l], we developed an adapted impulse sensor to nonintrusively investi- 
gate solitary wave propagation in a linear chain of identical elastic beads. We explored 
the problem of solitary wave reflection by changing the elastic properties of the wall and 
showed that the solitary wave detected at the wall differs from the actual solitary wave 
propagating through the chain. Our measurements significantly improve upon previous ex- 
perimental studies P, ^ and allows excellent agreement with our numerical simulations and 
Nesterenko's analytical theory 

The physical behavior of solitary waves in bead chains can be described as follows. Under 
elastic deformation, the energy stored at the contact between two elastic bodies submitted 
to an axial compression corresponds to the Hertz potential 0] Uh = (2/5)/t<5 5 / 2 , where 5 is 
the overlap deformation between bodies, kT 1 = (9 + 6»')(i? _1 + R'~ l ) 1/2 , 9 = 3(1 - z/ 2 )/(4F), 
and R and R' are radii of curvature at the contact. Y and v are Young's Modulus and 
Poisson's ratio, respectively. Since the force felt at the interface is the derivative of the 
potential with respect to 5, (F H = dsUn = n5 3 ^ 2 ), the dynamics of the chain of beads is 
described by the following system of iV coupled nonlinear equations, 



mdf t u n = k 



i \3/2 / n3/2 

K Un-i-U n )+ -{U n -U n+ i) + 



(1) 



where m is the mass, u n is the position of the center of mass of bead n, u n = 2nR at rest, 
and the label + on the brackets indicates that the Hertz force is zero when the beads lose 
contact. Under the long- wavelength approximation A ^> R (where A is the characteristic 
wavelength of the perturbation), the continuum limit of Eq. ^can be obtained by replacing 
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the discrete function u n ±i(t) by the Taylor expansion of the continuous function u(x±2R, t). 
Keeping terms of up to the fourth order spatial derivatives, Eq. [T] leads to the equation for 
the strain if) = —d x u > 0, 

Bfcl> * c 2 d 2 xx [^ 2 + (2/5)R 2 ^d 2 xx (^% (2) 

where c = (2R) 5 ^ 4 (K,/m) 1 ^ 2 jjj. Looking for progressive waves with speed v, in the 
form = x — vt), Eq. 121 admits an exact periodic solution in the form ip = (5/4) 2 x 
(f/c) 4 cos 4 [^/(i?vT0)] A Q, S, 0]. Although this solution only satisfies the truncated 
Eq. [TJ there is quantitative analysis on how well one hump (—ir/2 < £/(Ry/lQ) < ir/2) 
of this periodic function represents a soliton solution 0, Q] . Approximating the spatial 
derivative, the strain in the chain reads ip ~ 5/(2R), and the force felt at beads contact, 
F ~ K(2Rip) 3 / 2 , and v become, 



F m cos 6 



x — vt 



RVIO 
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In our experiment, we consider the chain of 21 identical beads of mass m, located on a 
Plexiglas linear track as shown on top of Fig. ^ A piezoelectric dynamic impulse sensor 
(PCB 208A11 with sensitivity 112.40 mV/N) located at the end of the chain provides the 
force at the rigid end. This sensor has a flat cap made of the same material as the beads. 
Beads are Tsubaki high carbon chrome hardened steel roll bearing (norm JIS SUJ2 equivalent 
to AISI 52100). The radius of the beads is R = 13 mm (tolerance is ±125 /im on diameter), 
and the density is p = 7780 kg/m 3 . The Young's modulus is Y = 203 ± 4 GPa [l^l, and the 
Poisson ratio is assumed to be v = 0.3; our beads have thus a k = 12 N/um 3//2 . Moreover, 

). Assuming that 

the contact surface is a disk of area A = n(9RF) 2 ^ 3 [15], the corresponding maximum 
compression force is roughly Fy ~ 470 N, which corresponds to an overlap 5y — 11 /jhi. 
Forces inside the chain are monitored by a flat dynamic impulse sensor (PCB 200B02 with 
sensitivity 11.24 mV/N) that is inserted inside one of the beads, cut in two parts. The total 
mass of the bead sensor system has been compensated to match the mass of an original bead. 
This system allows achieving non intrusive force measurement by preserving both contact 
and inertial properties of the bead-sensor system. The stiffness of the sensor k s = 1.9 kN//zm 
being greater than the stiffness of the Hertzian contact (k s 3> ku oc k5 1 ^ 2 ), means the 
coupling between the chain and the sensor is consequently negligible. To relate the force F s 



registered by the sensor with the actual force at the beads contact, we write the Newton's 
law for both masses, respectively, located in front (+) and in the back (— ) of the sensor. 
Thus, F± = F s ±m±dt t x±, with F s = k s (x + — This set of equations can be summarized 

as 

d 2 tt F s + uj 2 F s = ujI [(1 - 0)F+ + 0F-] , (4) 

where we have introduced the resonant angular frequency of the system u = \k s {m^} + 
m I 1 )] 1 / 2 , and the mass ratio (3 = m + /(m + + m_). Experimentally /3 = 0.11 and the resonant 
frequency, /o = ujq/(2tt) ~ 85 kHz, indicates that safe measurements can be obtained for 
signal whose period is greater than tq = 1/fo ~ 12 /is. However, a relation between F±(t) 
is needed to invert Eq. H]and then determine the force F + (t) or F_(t) from the force F s (t). 
Assuming that the pulse travels at a velocity v, this relation reads F-(t) = F + (t + t ), where 
to = (x+ — £-)/ v. An estimate of the velocity v is obtained from the time of flight of the pulse 
and the deconvolution of Eq. 0]by means of Fast Fourier Transform, then provides the actual 
force F + (t) felt exactly at the interface between two beads. Notice that the improvement 
introduced here represents a correction of the order of (3, i.e. about 10%. Signals from 
sensors are amplified by a conditioner (PCB 482A16), recorded by a two channels numeric 
oscilloscope ( Tektronix TDS340), and transferred to a computer. The acquisition is triggered 
by the contact between the small impacting bead and the chain; both being in contact with 
soft wires they cause the discharge of a capacitor in a resistor (1/RC ~ 1 /is). This circuit 
allows high repeatability, e.g. for time of flight measurements. In Fig. a solitary wave 
propagates along the chain of beads. The central peak corresponds to the impulse detected 
at the end, whereas the two peaks on the sides are the incident and reflected waves measured 
inside the chain. Notice that the central peak is much higher and broader than the actual 
solitary wave propagating along the chain, thus no quantitative information can be extracted 
from it without a detailed description of the interaction between the solitary wave and the 
wall. In order to characterize solitary waves, we look both for velocity and duration of 
incident pulses recorded at one contact far from the wall. According to Eq. El we map 
experiments to F(t) = F m cos 6 [(t — to)/r], to obtain the amplitude F m , the duration 2r, and 
the time of flight to of a pulse. To provide more accurate data for the velocity, we perform 
flight time measurements for different positions of the active bead. In addition, for every 
experimental configuration we record three sets of data and check repeatability, and the 
whole experiment is repeated three times. According to Eq. El we first look for the best fit 
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in a least squares sense for the experimental velocity of the pulse, in the form v = CF^ 6 , 
and we find an experimental value C e = 203.6 ±8.9 in standard units. This value agrees with 
the theoretical prediction C t = 198.9, derived from Eq. El within an error less than 3%. The 
fit is plotted in straight line in Fig. |2K- For sake of comparison, we also plot (the straight line 
on Fig. Et>) the duration 2r = 2RyW/v, also obtained from Eq. El The velocity is thus in a 
satisfactory agreement with the theoretical prediction, which also appears at first glance to 
predict in a good manner the duration of the pulse. However, energy dissipation is expected 
to produce a broader solitary wave. Dissipation is characterized by the restitution coefficient 
(see Fig. defined as e = (U n+ i/U n ) l/2 = (F n+1 /F n ) 5/6 (U n is the Hertz potential, i.e., 
the work done by the Hertz force F n at the contact n). Here we consider two mechanisms 
responsible for the dissipation; internal viscoelasticity and solid friction of beads submitted 
to their weight mg (g is the gravity), on the track. A third mechanism, the solid friction 
between beads due to thwarted rotations may also be taken into account. However, 
the contribution of a friction force of the form F* = h*kS 3 ^ 2 into Eq. ^ reduces simply 
to considering an equivalent nonlinear stiffness k* = (1 + Viscoelastic dissipation is 



included by using the simplest approximation [19J, |20j for which the dissipative force at the 
contact of two beads reads, F v — rindtiS 3 / 2 ), where 77 includes unknown coefficients due to 
internal friction of the material |l5l. l20j . Solid friction is taken into account by considering 
a frictional force F s = \img [l^. The potential energy difference (U n — U n+ x) being equal 
to the work done by both previous dissipative forces allows us to estimate the restitution 
coefficient to be force dependent, e = (U n+ i/U n ) l l 2 ~ 1 — C v F l l & — C s /F. Simple calculations 
provide the relation of 77 and /i with the new constants C v and C s as, 77 fa 2\/lORC v /5C 
and /i ~ AC s /5mg respectively. Experimentally, we determine that C v = 1.9 x 1CT 2 and 
C s = 1.7 x 10 _1 in standard units, see Fig. Efc. Then, r/ e w 1.8 fis and fi e ~ 0.19. 

Numerical simulations based on a Velo city- Ver let algorithm allow to explore the main 
features of solitary waves by solving Eq. directly. We first run numerical calculation 
without dissipation, plotted in dashed lines on Fig. |2K andEb- Looking for least square fit 
for the velocity, as previously done, we find C n = 201.5 ± 0.1. Compared to the theoretical 
value Ct, simulations improve the agreement with experiments (relative error on velocity is 
about 1%), but a noteworthy disagreement is now observed for the duration of the pulse (see 
Fig. Eh), which is about 10% lower than experimental values. This lag is consistent with 
the presence of a weak dissipation. At this stage, we only consider the effect of viscoelastic 
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dissipation in numerical simulations. We thus adjust the coefficients, and for rj n = 2 /is and 
H = 0, a good agreement can be obtained both for the velocity and the duration, in the range 
of amplitude where viscoelastic dissipation dominates over solid friction (F m > 50 N). Notice 
that r] n differs from the experimental value rj e only by 20%. Since solid friction has not yet 
been included in simulations, the experimental pulse is still broader than in simulations at 
low force amplitude (F m < 20 N) where this mechanism dominates. 

We now check how simulations reproduce the features of the reflection process. Figure 
shows the corresponding numerical simulations for the incident and the reflected solitary 
wave as well as the force registered at the wall. Although simulations include only viscous 
dissipation, rj n = 2 /is, the agreement between Fig. [Tk. and Fig. [TJd is very good. Notice that 
momentum is conserved, i.e., the area of the central peak in Fig.^i is equal to the area of the 
incident plus the reflected solitary wave. Figure^ presents the corresponding calculations 
of the time evolution of the potential and kinetic energy when a solitary wave interacts with 
the wall sensor. The solitary wave is initiated at t = by a purely kinetic impact. At 
t — 1 ms the pulse reaches the rigid sensor and the energy is stored into potential. The 
pulse is then reflected and propagates backward to the free end until leading to ejection of 
beads after t — 2 ms. 

We further investigated the solitary wave reflection by varying the mechanical properties 
of the flat part of sensor in contact with the last bead. This is done by locating polished 
disks of 1 mm thickness and 5 mm diameter of different known materials on the active part 
of the sensor. These samples are made of plexiglass, Mg, Cu, Si, Fe, and W. For materials 
softer than the beads, unexpected features arise. For instance, in Fig. Eb, the experimental 
force on the wall exhibits a well defined secondary peak. The break of symmetry implied 
by the change of elastic properties leads to the generation of a so-called secondary solitary 



121 ]. Dissipationless 



wave in the reflected impulse predicted recently via simulations in [ill , 
numerical simulations in Fig. Efc reproduce well the experimental finding of Fig. without 
adjustable parameter. Better agreement can be achieved but it requires the knowledge of 
the mechanism dominating dissipation of the samples. FigEK is the ratio of the maximum 
force measured at the wall and the respective maximum force of the incident solitary wave. 
Despite the peculiar form of the force, the ratio of maximum forces follows a well defined 
law that is characteristic of the kinetic to potential energy conversion at the wall. This 
interesting feature should prove valuable to determine the Young modulus of materials of 
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unknown nature, when the sample size is a practical limitation. 

To understand the underlying physics of solitary wave reflection, we focus on the kinetic- 
potential energy conversion when a solitary wave interacts with a rigid wall. As shown on 
Fig. ^ and P, Q, when solitary wave propagates freely in the chain, the kinetic energy 
Kchain is about 56% and the potential energy U c h a in is about 44% of the total energy (for 
a rough estimation we assume K c h a in — chain)- However, when a solitary wave reaches the 
end of the chain, the potential energy stored at the sensor-bead contact equals the total 
energy carried by the solitary wave. The kinetic energy is thus transformed into potential 
at the contact. Then, U^f ~ 2U c h a in- On the other hand, the solitary wave extends on a 
few beads, and the potential energy stored in the chain is roughly supported by the most 
compressed contact (U c h ain ~ U^d)- ^ Anally becomes, 

TTmax Tpmax / V \ — 2/5 

U end _.r, . r end p6/5 I , , 1 bead \ ,_x 

Tmax { ' v / v / 



TTmax Tpmax \ ~y 

u bead 1 bead v 1 end 

which is a function of Young modulus of beads and the sensor plane. In Fig. EK, we com- 
pare experiments, numerical simulation, and the above estimate. Within the error bars, a 
satisfactory agreement is obtained. 

In conclusion, we have developed a non intrusive reliable method to investigate solitary 
wave propagation and solitary wave reflection at walls. Our measurements in conjunction 
with our numerical simulations provide a powerful tool to accurately investigate a variety of 
related problems such as the main features of solitary waves reaching impedance mismatch, 
the generation of the recently predicted secondary solitary waves at the boundaries, and the 
solitary wave interactions, among others. 
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FIG. 1: (top) Schematic view of experimental setup, (a) Experiments: middle peak indicates the 
force signal at the end of the chain, whereas lateral peaks are the incident and reflected solitary 
wave. The solid line represents the force at a single contact extrapolated from Eq. |1J (b) and (c) 
The numerical simulations of the contact forces and energy, respectively, for r\ n = 2//s. 
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FIG. 2: (a) Velocity v and (b) duration 2r of the solitary wave, measured inside the chain, vs 
force amplitude. Theoretical predictions from Eq. |2J and numerical simulations are contrasted to 
experimental data, (c) Restitution coefficient vs force amplitude. 
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FIG. 3: (a) Ratio of force amplitudes at the end of the chain and of the incident impulse vs Young's 
modulus of the sample placed on the rigid sensor. Inset (b) force measurements when a solitary 
wave collides on the softer sample (Y = 5 GPa), and (c) corresponding simulation reproducing all 
the experimental features for r] n = 0. Dashed lines indicate forces at the end of the chain. The last 
peaks on the right represent the secondary solitary waves. 
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